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ABSTRACT 

We quantify the effects of magnetic fields, cosmic rays and gas pressure on the rotational velocity of H I gas in 
the Milky Way, at galactic distances between R Q and 2Rq. The magnetic field is modelled by two components; 
a mainly azimuthal magnetic component and a small-scale tangled field. We construct a range of plausible 
axisymmetric models consistent with the strength of the total magnetic field as inferred from radio synchrotron 
data. In a realistic Galactic disk, the pressure by turbulent motions, cosmic rays and the tangled turbulent field 
provide radial support to the disk. Large-scale (ordered) magnetic fields may or may not provide support to the 
disk, depending on the local radial gradient of the azimuthal field. We show that for observationally constrained 
models, magnetic forces cannot appreciably alter the tangential velocity of H I gas within a galactic distance of 
2R G . 

Subject headings: galaxies: halos — galaxies: kinematics and dynamics — galaxies: magnetic fields — galax- 
ies: spiral — dark matter 



1. INTRODUCTION 

The interstellar medium in galaxies contains three basic 
constituents: ordinary matter, cosmic rays and magnetic 
fields. Studies of the vertical distribution of gas and syn- 
chrotron emission in the solar neighbourhood show that cos- 
mic rays and magnetic fields influence the spatial distribu- 
tion of gas providing efficient support against the gravitational 
force (e.g., Ferriere 2001; Cox 2005). In the radial direction, 
gradients in the pressure may produce a difference between 
the rotational velocity of the gas and the real gravitational 
circular velocity v c , defined as = Rd$ / dR. Here R is the 
galactocentric radius and $ the gravitational potential. The 
asymmetric drift, defined as - measures this difference. 
In a gaseous disk in equilibrium, the asymmetric drift is a 
consequence of the support by thermal, turbulent and mag- 
netic pressures as well as the pressure due to cosmic rays (e.g., 
Parker 1966; Spitzer 1978). In galaxies with circular veloci- 
ties > 50 km s" 1 , the asymmetric drift corrections to derive 
the real gravitational circular velocity from the observed ro- 
tational velocity are not applied because they are small com- 
pared to uncertainties due to inclination, warps, non-circular 
motions, etc (e.g., de Blok & Bosma 2002). Only for low- 
mass galaxies with < 50 km s" 1 , corrections for the asym- 
metric drift must be taken into account (e.g. Dalcanton & Stilp 
2010). 

In this approach, magnetic effects on the gas are modelled 
as a pressure term in the asymmetric drift. However, gas can 
also experience an additional force due to the magnetic ten- 
sion of a large-scale magnetic field (e.g., Sanchez-Salcedo 
1997a; Ruiz-Granados et al. 2010, 2012; Jalocha et al. 2012). 
The magnetic tension by (ordered) large-scale magnetic fields 
may provide an inward or an outward force, depending on the 
local radial derivative of the azimuthal field (Sanchez-Salcedo 
& Reyes-Ruiz 2004). The energy density in the regular mag- 
netic field might be more important than has often been be- 
lieved (e.g., Beck 2007). For instance, Fletcher et al. (2004) 
find that the regular magnetic field in M3 1 is almost perfectly 
axisymmetric and essentially constant (~ 5/iG) between 8 and 



14 kpc. 

Recently, Ruiz-Granados et al. (2010, 2012) claim that 
large-scale magnetic fields can provide enough radial confine- 
ment of the gas to explain the rising-up in the H I rotation 
curve detected in some galaxies. Moreover, they argue that 
the shape of the H I rotation curves of M3 1 and the Milky Way 
are fitted better if the contribution of the large-scale (mainly 
azimuthal) magnetic field is included. It remains unclear how 
these findings are compatible with rigorous upper limits based 
on the Virial Theorem arguments indicating that magnetic 
fields can hardly speed up H I disks by more than 20 km s -1 
in the outermost point of H I detection (Sanchez-Salcedo & 
Reyes-Ruiz 2004). Since the relative importance of turbulent, 
magnetic and cosmic ray pressures is comparable, it is cer- 
tainly not clear what is the role of the pressure by the cosmic 
rays and the tangled component of the magnetic field in pro- 
viding support to the disk. In this paper, we combine different 
observations to explore on whether and how magnetic fields 
and cosmic rays can alter the gas rotation curve in the Milky 
Way. 

The paper is structured as follows. In Section [2] we de- 
scribe the formalism, our simplifying assumptions and the 
basic equations. In Section [3j we highlight the reasons why 
the Milky Way is an excellent target to carry out this analysis 
and present a range of plausible magnetic models compati- 
ble with observations. In Section|4]we calculate the expected 
differences between the observed rotation curve and the true 
gravitational circular velocity for these models. Conclusions 
are given in Section [5] 

2. ASSUMPTIONS AND GOVERNING EQUATIONS 

We consider a magnetized disk of slightly ionised gas with 
the axis along the z-direction that is described well in the 
ideal magnetohydrodynamic limit (MHD). In galactic disks, 
the cosmic ray population forms a light fluid with significant 
pressure, which is coupled via magnetic fields to the thermal 
interstellar components. Hence, the pressure by cosmic rays 
may help to support thermal gas and will be also included. 

Following previous works, we will assume that the disk 
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FIG. 1 . — Midplane H I volume density of the Milky Way taken from KD 
(solid line) and NS (dashed line). The KD curve corresponds to an exponen- 
tial fit to the data. 



is axisymmetric over time (Sanchez-Salcedo 1997a; Jalocha 
et al. 2012; Ruiz-Granados et al. 2010, 2012). Non- 
axisymmetric configurations are more difficult to treat be- 
cause they generate magnetic density waves (e.g., Lou & Fan 
1998). Axisymmetry is the simplest assumption to quantify 
the overall effect of magnetic fields on the azimuthally aver- 
aged tangential velocity in an equilibrium configuration. 

It is assumed that the magnetic field can be decomposed 
into an average part B(R) varying only on the large scale and 
a small-scale isotropic random field b, so that (b) = 0. We 
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will refer to (b 2 ) as the strength of the random (or tur- 
bulent) magnetic field. At scales larger than the coherence 
length of the small-scale magnetic field, it is useful to define 
the strength of the total magnetic field as B 2 0t = B 2 + (p 2 ). 
In the equilibrium configuration, we assume that the_ regu- 
lar magnetic field consists of a planar magnetic field B(R) = 
(Br,B0,O) (in cylindrical coordinates), with Br <C B$, We fur- 
ther assume that the radial velocity of the gas, Vr, is much 
smaller than v^, and thus it can be ignored; thus the veloc- 
ity in the disk is v = (Q,v^,Q) in cylindrical coordinates 1 . In 
principle, each component of the interstellar gas can rotate at 
different velocity. Since we are only interested in the rotation 
curve of neutral atomic gas, we will consider the dynamics 
of this component and ignore the presence of molecular hy- 
drogen gas. In the Milky Way, this is a good approximation 
especially at R > 10 kpc because it is at these galactic dis- 
tances where the neutral atomic hydrogen is dominant in the 
mass budget of the interstellar gas. 

Because of the symmetry about z = 0, we take that all the 
derivatives with respect to z are negligible near the midplane 
of the disk. Under these circumstances, the radial component 
of the motion equation of the gas at z = reads 



v l = v2 c +v p +v L g , 



(i) 



where v^ ag is the contribution of the regular (azimuthal) mag- 
netic field: 



1 d(R 2 Bp 
8nRp dR 




(2) 



1 This is an approximation because the magnetic field creates a torque, 
unless Br = 0, leading to a radial inflow of gas (Sanchez-Salcedo 1997b). 



and vj, is the contribution by pressure gradients, 

RdP T 



..2 _ 



p dR 



(3) 



where p is the gas volume density at the midplane and Pt(R) = 
Pg+Pb+PcR is the total gas pressure consisting of the gas ki- 
netic pressure (thermal plus turbulent), the magnetic pressure 
Pb, arising from the random magnetic field component plus 
also the pressure by cosmic rays Pqr (the pressure by radia- 
tion will be ignored). More specifically, the kinetic pressure 
P g is given by 



P g = P<J, 



(4) 



where a is the H I line width in the radial direction, which is 
approximately constant or slightly decreasing with R in the 
outer parts of the H I disks, typically a ~ 6 - 8 km s" 1 (e.g., 
Dib et al. 2006, and references therein; see Blitz & Spergel 
1991 and Burton 1992, for our Galaxy). The spatially aver- 
aged magnetic pressure by the turbulent field is taken as 



(5) 



Finally, the cosmic-ray pressure is expected to be proportional 
to magnetic pressure: 



P CR - 



Mg|t 
8tt ' 



(6) 



where B loX is the strength of the total magnetic field (that is, 
B 2 ol = B 2 + {b 2 )) and p is a constant of the order of 1. This 
is justified by minimum-energy-type arguments (e.g., Beck et 
al. 1996). 

We want to stress that vj, and v^ ag are not necessarily posi- 
tive quantities. For instance, an unmagnetized isothermal disk 
with a radially decreasing density has vj < 0, which signifies 
that it provides pressure support to the disk (i.e. < v c ) be- 
cause it produces a force pointing outward. From Equation 
(0, it is simple to see that the magnetic tension imposes an 
inward force, i.e. v^ ao > 0, provided that the azimuthal mag- 
netic field decays radially not faster than l/R. Note that in the 
axisymmetric case with B z = 0, the divergence-free condition 
implies Br oc l/R. Therefore, if the magnetic pitch angle is 
constant with R, we infer B^ oc 1 /R, implying that v^ ag = 0. 
Consequently, a radial decay of B^ slower than l/R requires 
a pitch angle decreasing with R. 

To study the distribution of mass of a certain galaxy, we 
need v c , that is the circular speed of a test particle, but what 
we observe is the azimuthal velocity of the gas v^. For what 
values of v|, + v^ ag is the correction to the rotation curve sig- 
nificant? As guide numbers, if we observe that the gas rotates 
at a given galactocentric radius with a tangential velocity of 
= 230 km s" 1 and vj, + v 2 rag ~ 9000 km 2 s" 2 , then the grav- 
itational circular velocity is v c ~ 210 km s -1 . Hence, val- 
ues of 9000 km 2 s" 2 produce a boost of 20 km s" 1 . On the 
other hand, if vj, + v 2 was -9000 km 2 s" 2 , then v c = 249 km 
s . For a low-mass galaxy with = 120 km s , a value of 
: 4400 km 2 s" 2 implies that v c = 100 km s" 1 . It is 



"P T ''mag 

v p + v mag tnat we want to calculate in the Milky Way. 



3. THE MILKY WAY AS A TEST CASE: MODELS 

The determination of the distribution of H I volume den- 
sity and the magnetic structure of our Galaxy has been im- 
proved significantly over the last two decades. Therefore, our 
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FIG. 2. — Variation of the total (Btot; dot-dashed lines), regular (B^; solid lines) and turbulent ((fo 2 ) , dashed lines) magnetic field strengths for the different 
models. In this plot, we have used Rq = 8.5 kpc and the reference KD model. The symbols show the strength of the total magnetic field from the deconvolved 
surface brightness of synchrotron emission at 408 MHz (see text). 

cm -3 and Rh = 3.15 kpc. For comparison, we also plot the ex- 
ponential fit as derived by KD. in Figure Q] We see that the 
major discrepancies between NS and KD occur inside 1 .5Rq. 
KD derived an Hi plateau in surface density of 10M©pc" 2 at 
the inner Galaxy, which fits better what we know for external 
galaxies than the saturation value derived in NS, of 2M Q pc~ 2 . 
Therefore, we will use KD as our reference H I gas model but 
also refer to NS to explore the effect of systematic uncertain- 
ties in the the derivation of the midplane H I density. In order 
to include 28% of helium and 1.5% heavier elements, we will 
convert n# in mass density by the relation p=l Am p riH, where 
m p is the proton rest mass (e.g., Ferriere 2001; Cox 2005). 

The magnetic field of the Galaxy has been studied through 
synchrotron emission, Faraday rotation, optical polarization 
and Zeeman splitting. The strength of the total magnetic 
field averaged in azimuth, Z? tot , was obtained from the surface 
brightness of synchrotron emission at 408 MHz. In the radial 
interval between 3.5 kpc and 17 kpc (using Rq = 8.5 kpc), 
and assuming energy equipartition between magnetic fields 
and cosmic rays, the strength of the total magnetic field in the 
disk can be fitted by 



Galaxy is a natural laboratory to quantify the effects of mag- 
netic fields in the gas dynamics (see also Vallee 1994). In or- 
der to estimate v 2 P and v 2 ™ we need to know the azimuthally 
averaged radial distribution of the H I volume density at the 
midplane, and the radial profile of both and {b 2 ). 

The azimuthally averaged H I volume density at the mid- 
plane has been derived by Nakanishi & Sofue (2003, here- 
after NS) using the Leiden/Dwingeloo survey, the Parkes 
survey and the NRAO survey. More recently, Kalberla & 
Dedes (2008, hereafter KD) inferred the average H I den- 
sity at the midplane excluding extra-planar gas, using the 
Leiden/ Argentine/Bonn (LAB) H I line survey, which com- 
bines the southern sky survey of the Instituto Argentino de 
Radioastronomfa (IAR) with an improved version of the Lei- 
den/Dwingeloo survey. NS and KD used different assump- 
tions to derive the H I structure of the Milky Way. NS chose 
as Galactic constants Rq = 8 kpc and v© =217 km s" 1 , and 
adopted a slightly declining rotation curve to convert the ob- 
served brightness temperature distribution to density, whereas 
KD used/? Q = 8.5 kpc and v Q = 220 km s -1 and an almost flat 
rotation curve. In addition, NS assumed cylindrical rotation 
along z, whereas KD adopted a "lagging" halo. A discussion 
about the impact of the different assumptions on the H I dis- 
tribution can be found in Kalberla et al. (2007). 

In Figure Q] we plot the midplane H I density distributions 
versus R/Rq, as derived in NS. We see that the density de- 
cays exponentially beyond 1.5^ Q . On the other hand, KD 
found that, for 7 < R < 35 kpc, the midplane H I density can 
be approximated by n# = « Q exp[-(/?-/? Q ) /Rh] with n = 0.9 



•Btot.fit = B tot: Q exp 



R-R, 



R n 



(7) 



where Z? to t.© is the total magnetic field near the Sun, which is 
about 6 ± 2 /J3, and R B = 12 kpc (Beck et al. 1996; Strong et 
al. 2000; Beck 2001; Ferriere 2001; see also Jansson & Farrar 
2012 using WMAP7 22 GHz data). Since there is no reliable 
observational measurement of the magnetic field strength be- 
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FIG. 3. — Magnetic contribution to the rotation curve due to the azimuthal 
magnetic field, v raag , as a function of radius, for models A, Bl, B2 and C, 
using the KD profile (upper panel) and the NS profile (lower panel). The gas 
rotates at a speed of : 
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yond 2Rq, we will restrict our analysis to R < 2Rq, the same 
interval as studied in Ruiz-Granados et al. (2012). 

In order to estimate v^ ag and Pb we need to separate the 
ordered magnetic field and the turbulent magnetic field com- 
ponents. Starlight and synchrotron polarization data suggest 
that the local ratio between the regular and the total magnetic 
fields is 0.6-0.7 (e.g., Berkhuijsen 1971; Brouw and Spoel- 
stra 1976; Heiles 1996; Beck 2001). This implies that the 
regular magnetic field is 4 ± 1 /iG at the Solar radius. On 
the other hand, from Faraday rotation of pulsars and radio 
sources, Han et al. (2006) derived a regular field strength of 
2.1 ±0.3 /iG at the Sun position. Possible explanations for the 
difference between the equipartition estimate and the value in- 
ferred from pulsar data were discussed in Heiles (1996) and 
Beck et al. (2003). Since our aim is to place upper limits on 
the magnetic effects, we will take the value derived from po- 
larization measurements, Bq = B^Rq) ~ 4 /iG, as a generous 
value. 

The radial profile of B^ is not well constrained by obser- 
vations. In the inner Galaxy (3 kpc < R < Rq), the ordered 
magnetic field gets stronger at smaller Galactocentric radius, 
probably as R~ l or RT 2 (Heiles 1996). Han et al. (2006) used 
an exponential function to fit the ordered magnetic field and 
found a scale radius of 8 .5 ± 4.7 kpc in the radial interval be- 
tween 3 kpc and Rq. In the outer Galaxy (R > Rq), there 
exists no quantitative estimate of its exact R dependence, ex- 
cept that Z?0 < B lot . As already said, if the magnetic pitch 
angle is assumed to be constant with R, then B^ oc R~ l . At 
R < 2R Q , this radial decay is consistent with WMAP7 22 GHz 
data (Jansson & Farrar 2012). 



In order to illustrate how the results depend on the assump- 
tions, we will explore four different representative magnetic 
configurations (see Figure |2). In model A, we will assume 
that declines exponentially with R, in the outer Galaxy 
(R > Rq), with the same scalelength as B tot : 

R-Rq 



fi (fl) = fl exp 



Rr 



(8) 



with B Q = 4 /iG and Rg = 12 kpc. As a consequence, the ra- 
tio of the regular magnetic field to the total magnetic field, r\, 
is constant with radius for R > Rq, having a value of ~ 0.6. 
This is a well motivated possibility because constant values 
for 7/ with galactocentric distance have been derived in ex- 
ternal galaxies. For instance, in the case of M31, Fletcher 
et al. (2004) derived r\ ~ 0.7 in the radial range 8 to 14 kpc. 
A rather constant value of r\ within R < 6 kpc was found by 
Beck (2007) for the galaxy NGC 6946. In M33, Tabatabaei et 
al. (2008) inferred a value of 77 ~ 0.45 independent of radius 
within R < 7 kpc. 

In a second type of magnetic profiles (labeled as models B 1 
and B2), we will adopt the same dependence for the azimuthal 
field of our Galaxy as in Ruiz-Granados et al. (2012): 



B^{R) = 



(R, + Rq)Bq 
R, + R ' 



(9) 



where Bq = 4 /iG and R/ is, in principle, a free parameter. To 
facilitate comparison with previous work, we will take Ri = 14 
kpc. Models Bl and B2 have the same regular magnetic field 
as given in Eq. (O but differ in the random magnetic compo- 
nent. In model Bl, we will assume that rj is constant with ra- 
dius, and has a value of 0.66, rather similar to model A. Thus, 
(b 2 ^ = (?/~ 2 - 1)^0 ~ 1-3B^- In model B2, the mean square 

turbulent field (Z? 2 ) is obtained as the difference between the 
total magnetic field as inferred from synchrotron emission and 
the regular magnetic field, that is (b 2 } = B^^-B^. where 
#tot,fit is given in Eq. ^} with B to t,© = 6.9 /iG. 

Finally, we consider a fourth model (labeled as model C) 
in which we assume that there is equipartition between the 
magnetic pressure in the random field and the dynamical pres- 
sure, that is P/, = P g at any radius in the range Rq < R < 2Rq. 
This is expected in turbulent disks where turbulent motions 
in the gas tangle the magnetic field. The equipartition con- 
dition determines (b 2 ^ as a function of radius. Once (b 2 } is 
derived, the coherent magnetic field B$ is then obtained as 
B~l=Bl fit -(b 2 ). 

Figure |2] shows the radial profiles for both the strength of 
the azimuthal large-scale magnetic field and the strength of 

1 11 

the small-scale random field (b 2 ^ , for the different models. 
We have assumed that Rq = 8.5 kpc. Note that the magnetic 
profiles in model C depend on the adopted midplane H I den- 
sity; to easy discussion we show the magnetic profiles for our 
reference KD density profile. For comparison we also plot 
the total magnetic field as derived from synchrotron emission 
(Beck et al. 1996; Beck 2001; Ferriere 2001). By construc- 
tion, models A, B2 and C fit very well the strength of the total 
magnetic field at R > 7 kpc. Bearing in mind that the error 
in the determination of B tot from synchrotron emission data 
is 30%, we can see that models Bl are compatible with it at 
R > 7 kpc. 

In model A, the azimuthal magnetic field varies between 
4/iG at Rq to 2.0/iG at 2Rq, and has 77 ~ 0.6 constant with 
radius. In model Bl, B^ varies between 4/iG and 2.9 /iG, and 
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FIG. 4. — Difference between and v 2 c for models A, Bl, B2 and C, and different combinations of a and fj,. For the midplane density we have used the 
reference KD profile. A negative (positive) value means that the gas rotates slower (faster) than the gravitational circular velocity. Cases with cr = corresponds 
to ignoring the kinetic pressure of the gas, whereas models with fj, = assume that the radial force by the pressure of cosmic rays is null. 

simple to show that, in model A, v^ ag < at R > Rb- On the 



r] is also constant with radius 77 ~ 0.66. Model B2 has the 
same azimuthal magnetic field as model Bl but 77 increases 
radially from 0.58 at R Q to 0.9 at 2R Q . In model C plus the 
KD density profile, the azimuthal magnetic field varies be- 
tween 4.6 pG at Rq to 3.1 pG at 2Rq, and 77 increases from 
0.68 to 0.9. Finally, in model C plus the NS density profile, 
Bcf, varies between 6.3 pG at Rq to 3.4 /iG at 2Rq, and 77 > 0.9 
at any radius betweeen R Q and 2R Q . Thus, the magnetic field 
at 2Rq is dominated by the regular component in models B2 
andC. 

We should note here that the different radial profiles for 
are realistic for a finite radial interval, Rq < R < 2R Q , but 
they may not so realistic at large R. For instance, the mag- 
netic field strength in models Bl and B2 declines asymptoti- 
cally as R~ l at large R, which is unphysical because, accord- 
ing to Ampere's law, this requires a net electric flux through 
the equatorial plane of the galaxy (Sanchez-Salcedo & Reyes- 
Ruiz 2004). 

4. ESTIMATING THE CONTRIBUTIONS TO THE ROTATION CURVE 
4.1. Results 

We will start our discussion by considering the effect of the 
azimuthal magnetic field in the rotation curve, v^ ag . Figure 

[3] shows v^ ag as a function of R for the different models. We 

see that the shape of v^ a „ as a function of R depends critically 

on the adopted profile for B^. In models A and C, v 2 is 
positive at small galactocentric radii but turns out to negative 
values beyond a certain radius. Using Eqs. (0 and (|8), it is 



other hand, 



*mag i s positive at any radius in models Bl and 
B2. The sign of v^ a(J at 2Rq is model dependent and there is 
no a clear preference for a model with v 2 > over another 

with v^ ag < at 2Rq. Since, according to Equations (O and 
(0, the strength of the radial force (per unit of mass) by mag- 
netic effects and cosmic rays is proportional to p~ l , the exact 
values for v 2 and v\ depend on the H I gas model. If the 
NS profile is used, v^ ag at 2Rq ranges between -200 km 2 s~ 2 

to 400 km 2 s" 2 depending on the model, whereas it varies be- 
tween -100 km 2 s" 2 to 220 km 2 s~ 2 when the KD profile is 
used. Therefore, it has a minor effect on the rotational veloc- 
ity of the gas; the corresponding correction is ~ 0.5 - 1 km 
s" 1 at 2Rq. At this galactocentric radius, the correction by 
v 2 is comparable to the correction by the kinetic pressure of 
the gas. For instance, consider the H I gas model of KD. The 
pressure correction is 



RdP g R 2 



(10) 



p dR R H 

which is ~ -260 km 2 s" 2 at 2R Q (using Rh = 3.15 kpc and 
<T = 7 km s" 1 , see $2] and $3). In the next, we discuss and 
quantify the relative importance of v 2 P as compared to v 2 . 
Figures |4] and show the contributions to the tangential ve- 



locity of the gas, vi 



, when the magnetic pressure P/, is in- 



cluded, for different combinations of P g and P CR . Obviously, 
the curves with a = and p = correspond to P g = Pqr = 
[see Eqs. © and ©]. On the other hand, curves with a = 7 
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FIG. 5. — Same as Figure[4]t>ut for the NS gas profile. For consistency, we use Rq = 8 kpc here. 



km s" 1 and /i = 0, include the kinetic pressure of the gas and 
the magnetic forces (i.e., including both the azimuthal and the 
small-scale components), but not the pressure by cosmic rays. 
The case /i = 1 describes a situation in which the pressure by 
cosmic rays is in equipartition with the magnetic pressure. For 
instance, Ferriere (2001) quotes a midplane value of [i = 1.28 
in the vicinity of the Sun. In order to interpret correctly Fig- 
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FIG. 6. — Total gravitational circular velocity v c (solid line) together with 
the tangential velocity of the gas vj, after including magnetic fields, gas pres- 
sure and cosmic rays in model A (dotted lines). The upper dotted line was 
derived using the H I density profile from KD, whereas the lower dotted line 
was calculated using the profile derived in NS. Symbols indicate the observed 
rotational velocity of Hi gas taken from Ruiz-Granados et al. (2012). The 
contribution of the different mass components to the rotation curve are also 
shown: bulge (dot-dashed line), stellar disk (long dashed line), gas (short 
dashed line) and dark halo (triple dot-dashed line). We see that the difference 
between and v c is small as compared to the observational uncertainties. 



ures |4] and [5] remind that when 



is negative, it means 



that the MHD terms provide support to the disk and, hence, 
the measured tangential velocity lags the circular velocity of 
a test particle, i.e. < v c . 

In our models, the strength of the random field (b 2 ) 1 ^ 2 de- 
clines with radius and, therefore, the magnetic pressure by this 
small-scale magnetic field produces a force outwards, giving 
support to the disk (slower rotation). In models A, B2 and 
C, this outward force is able to compensate any pinching ef- 



fect by the azimuthal magnetic field. Thus v j 



< at any 



galactic radius > 7 kpc. 

As can be seen in Figure |4] and [5] the effect of magnetic 
fields in the rotation curve is less important in model Bl, even 
though it presents the highest B tot -values at 2R Q (Figure |2). 

_ 1/2. 

The reason is that the radial profiles of both and lb 2 ') 
in model B 1 are more shallow and, thus, the confining effect 
of Bfj, is partly balanced by the radial support of the random 
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FIG. 7. — Azimuthal magnetic field and the corresponding effect on the 
gas rotation velocity for a double piece-wise exponential magnetic field, with 
r\ = 0.6, a = 1 km s -1 and fj, = 1, all constant with radius. In all the models, 
the magnetic field declines exponentially with a scalelength of 12 kpc up to 
2Rq . Beyond 2Rq kpc, the scalelength of the magnetic field is 3 kpc (dotted 
line), 6 kpc (dot-dashed line) and 12 kpc (dashed line). The mass model is 
the same as in Figure|6] Symbols indicate the observed rotational velocity of 
Hi gas taken from Ruiz-Granados et al. (2012). Note that the y-axis ranges 
from 100 to 260 km s -1 . 

fields. In model Bl, the outward force by the radial gradient 
of Pj, is in balance with the inward pinching force created by 
Z?0 at a radius 

„2 



1- 



■Tf 



(ID 



For Ri = 14 kpc and r/ = 0.66, this balance occurs at R = 10.8 
kpc. In order to have Rb\ > 20 kpc in this kind of models, we 
need 77 > 0.76. In model B2, the equality between the gradient 
of Pb and the inward force by occurs in a more inner radius. 
Model B2 illustrates the role of the magnetic pressure by the 
random field; even if 77 = 0.9 at 2Rq, the contribution of the 



which is Rp dPb/dR, is three times 



random field to 

larger than v 2 ^. 

As already mentioned, the contribution of the kinetic pres- 
sure of the gas to the rotation velocity is comparable in magni- 
tude to the magnetic terms and, therefore, it must be included 
to estimate correctly the asymmetric drift v^-v 2 . The same 
holds true for the cosmic ray pressure. In our models, when 
all the terms are added, - v 2 < at any radius of inter- 
est. Hence the gas rotates at a speed less than v c . For the 
KD density profile with a = 7 km s" 1 and /x = 1, we infer that 
v p +v mae at ranges between -1000 to -1850 km 2 s~ 2 , and 



between -1500 to -2600 km 2 s" 2 if the NS density is used. 
Figure [6] shows the gravitational circular velocity for a typ- 
ical mass model of the Galaxy, together with the tangential 
velocity of the gas in model A, when all the correction terms 
in the asymmetric drift are included. We see that the effect 
is very small as compared to the intrinsic uncertainties in the 
determination of the H I rotational velocity; at 2Rq the gas is 
expected to rotate about < 4- 8 km s" 1 slower than the corre- 
sponding gravitational circular velocity v c at this radius. We 
find that for realistic models, the magnetic fields and cosmic 
rays cannot rise the rotation curve at galactocentric distances 
of ~ 2Rq, contrary to the suggestion by Ruiz-Granados et 
al. (2012). 

It is interesting to note that Figures |4] and [5] show that the 
correction to the rotation curve by magnetic fields and cosmic 
rays increases steeply with galactocentric radius. In partic- 
ular, if we extrapolate model A to larger galactocentric dis- 
tances and use the NS gas profile, we would find that the gas 
at 3Rq would rotate ~ 40 km s" 1 slower than a test particle on 
a circular orbit (see Figure|7). However, the magnetic profiles 
used in our models are based on synchrotron observations at 
R < 2Rq. Thus, there is no reason to assume that these pro- 
files are valid at any R. In order to explore how depends 
on the adopted magnetic profile, Figure|7]shows in models 
with r\ = 0.6 (constant with radius), a = 7 km s" 1 and /i = 1, 
where the azimuthal magnetic field is described by a dou- 
ble piece-wise exponential profile; Rb = 12 kpc at R < 2Rq 
kpc (as in model A) but having a steeper radial decline be- 
yond 2Rq. We see that when Rb — 6 kpc in the outer disk 
(R > 2Rq), the MHD terms produce a shift in the azimuthal 
velocity of ~ 15 km s" 1 (the KD density profile was used). 
On the other hand, if the outer magnetic field declines with 
a scalelength of 3 kpc, the radial support by cosmic rays and 
magnetic fields leads to v c — > 5 km s" 1 only between 2Rq 
and 2.6Rq. Beyond 2R Q , empirical determinations of the 
strength/topology of magnetic fields and the H I rotation curve 
are challenging and, hence, there is no means of testing the ef- 
fect of magnetic fields on the H I azimuthal velocity. Still, one 
has to consider the vertical confinement of the magnetic fields. 
Consider model A with a single exponential scalelength of 
Rb = 12 kpc. At 3Rq, the cosmic ray plus magnetic pressure 
in the midplane is ~ 0.2 x 10~ 12 dyn cm" 2 . Since the weight 
of neutral gas cannot account for this large pressure, an addi- 
tional coronal component should be invoked to provide verti- 
cal support. Following the same analysis as Cox (2005) did 
at the solar neighbourhood, we find that a component with 
density (0.003cm~ 3 )exp(-|z|/z ? ), with z g ~ 10- 12 kpc, and 
temperature of ~ 3 x 10 5 K (at 3Rq) could confine the cos- 
mic rays and magnetic field at 3Rq. This coronal layer would 
probably produce excessive X-ray emission and the total col- 
umn density of OVI would be orders of magnitude more than 
is observed looking out of the galactic plane (Cox 2005). It 
is more simple to assume that beyond the stellar truncation 
radius, stellar formation is almost nonexistent, so energy in 
cosmic rays and magnetic fields decays faster with R because 
the energy input into these components from stellar processes 
is likely to be less important (e.g., Oiling & Merrifield 2000). 

4.2. Other magnetic models 

We have shown that in all our models, the magnetic fields 
provide support to the disk when the magnetic pressure by the 
random field component is included (at least beyond a galac- 
tocentric distance of 10.8 kpc). If we had adopted = 2/jG 
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instead of ApG, the magnetic pressure by the random field 
would increase by a factor of 1 .6 (in order to account for the 
synchrotron radio emission), providing more radial support to 
the disk. Consider, for instance, model A with the KD profile 
for a = 7 km s" 1 and p = 1. At 2Rq, vp + v^^ would change 
from -1850 km" 2 s" 2 for B Q = 4pG, to -2250 km 2 s" 2 if a 
lower value Bq = 2/xG is adopted. 

In order for the gas to rotate faster than a test particle, we 
need vp + v 2 > 0. This is possible only in models in which 
the radial gradients in the pressure created by the random 
magnetic field and cosmic rays are small. To explore this pos- 
sibility, we have treated R/ in model B 1 as a free parameter 
and allowed to vary until vp+ v 2 ^ reaches a maximum at 2Rq. 
We found that the maximum occurs when the radial scale 
length is very large so that B^, Pb and Pqr are constant with 
galactic radius. In this configuration, B tot ~ 6 pG, constant in 
the range Rq < R < 2Rq, which is inconsistent with the ex- 
ponential radial decline observed in the synchrotron emission. 
In addition, a constant pressure of cosmic rays can hardly be 
maintained if the sources of cosmic rays are related to star- 
forming regions in the inner galaxy and they diffuse outwards 
lossing energy (e.g., Everett et al. 2010). Still, it is worthwhile 
to consider what happens in this unlikely situation. Since the 
pressure gradients by the turbulent magnetic field and cosmic 
rays are null, only the kinetic pressure and the v^-term play 
a role. Combining both contributions and assuming that p de- 
clines exponentially with R, we find 



v p + v„ 



B% (R-Rq 
exp ' 



47T/9 



© 



R 



H 



T al - (12) 



For Bq = 4 pG, Rh = 3.15 kpc and a = 7 km s we infer 
v p + v mag = 637 km 2 s~ 2 . The correction to the tangential ve- 
locity of the gas is 1 .5 km s" 1 if v c = 220 km s" 1 or 1 .8 km s" 1 
if v c = 175 km s" 1 . We conclude that, in axisymmetric mod- 
els, the observed tangential velocity is not expected to differ 
significantly from the true gravitational circular velocity at the 
interval R Q < R < 2R Q . 

4.3. Comparison with previous works 

Ruiz-Granados et al. (2012) claim that a significant im- 
provement of the fit to the rotation curve of the Milky Way is 
obtained when magnetic fields are considered. These authors 
use the same expression for the radial profile of the azimuthal 
magnetic field as that given in Eq. (0 and find the values of R/ 
that provide the best fit to the shape of the rotation curves of 
the Milky Way, ignoring any contribution from the turbulent 



component of the magnetic field, kinetic pressure or cosmic 
rays. They find Ri = 14.2^'j} kpc in a mass model with a 
pseudo-isothermal dark halo and R/ = 16.5 ±1.1 kpc if there 
is no dark matter at all within a sphere of radius 2Rq. 

The exact values for Bq, p Q , Rq and Rh in Ruiz-Granados 
et al. (2012) differ from those adopted in this paper, but only 
slightly. They used Bq = 3 pG, Rq = 8 kpc, Rh = 4 kpc, a 
column density of gas at the Sun position of lOM pc" 2 , and a 
constant vertical scale height of 0.2 kpc across the disk. Using 
these values and Ri = 14.2 kpc, we obtain that v 2 = 75 kirr 

s" 2 , at a galactic distance of 2Rq. For/?/ = 16.5 kpc, v 2 ^ = 90 

km 2 s~ 2 again at 2Rq. Even if we neglect the radial support by 
the cosmic-ray pressure and by the turbulent magnetic field, 
and only include the kinetic pressure of the gas, we infer vp + 
v 2 ^ ~ -100 km 2 s~ 2 at 2Rq, which is unable to provide the 
desired effect in the rotation curve. 

5. CONCLUSIONS 

How magnetic effects alter the overall rotation curve of gas 
in galaxies is a reoccurring theme in the literature. In a re- 
cent paper, Ruiz-Granados et al. (2012) claim that magnetic 
field forces provide the simplest way to explain the peculiar 
rising-up of the rotation curve in our Galaxy. Our Galaxy of- 
fers a unique opportunity for studying the three-dimensional 
distribution of neutral gas and magnetic fields in a detail un- 
obtainable in external galaxies. We have explored a range 
of plausible models to quantify the contribution to the radial 
support by kinetic gas pressure, magnetic fields, and cosmic 
rays. We restrict ourselves to the interval Rq < R < 2Rq, be- 
cause beyond 2Rq there is no determinations of the strength 
of the magnetic field. We have shown that, even adopting 
magnetic field configurations with a regular field of ~ 3pG 
at 2Rq, the rotation curve of our Galaxy is not appreciably 
altered by magnetic effects. Turbulent motions, cosmic rays 
and the random small-scale component of the galactic mag- 
netic fields act as pressure, giving support to the disk and, 
therefore, leading to a rotation a few km s" 1 slower than the 
gravitational circular speed. Given the large uncertainties in 
the rotation speed of the outer parts of the Galaxy, of ±25 km 
s" 1 , they can be safely ignored at least within R < 2Rq. 
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